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Introduction & Motivation 


=Most nonlinear optimization methods require an initial guess “close” to 
the optimal solution 


"Initial guess can be hard to find 


"There may be many local optima, and the initial guess dictates which 
one is found 


"General NLP “black box” solvers (IPOPT, SNOPT, ... ) can be slow 
= Many iterations 


= For large problems, most CPU time is spent in the optimizer 


Solution Overview 


=Start with very poor initial guess 


"Use multiple shooting with ~50-100 nodes 
=" Continuous thrust during propagation 
= Constrain defects to go to zero at matchpoints 


"Optimize states, controls, and endpoint locations by solving 
a series of quadratic sub-problems 


Multiple Shooting Formulation 


"Trajectory defined by a set of nodes with 
= Position 


= Velocity 
=" Mass 
=" Continuous control 


=Propagate backwards and forwards in time 


"Enforce continuity at matchpoints 
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Multiple Shooting Formulation 


"Fixed step size integrator is used, with 5 steps per node 


=Why fixed step size? 
= More consistent finite-differenced partial derivatives > faster 
convergence 


= Faster integration (don’t get stuck at a singularity with poor initial 
guess) 


= Better for parallelization (future work) 


"Runge-Kutta 78 numerical integration is used 


= Normally, use the 8" order truncation term to estimate the error in 
the 7" order step. Then choose the largest step size possible where 
the error remains within tolerance. 


= Here, we force a fixed step size, but use the truncation term to 
output the error estimate for use in mesh refinement 


=Mesh refinement: Add nodes where the 8" order 
truncation term for any of the integrator steps is > tolerance 
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Multiple Shooting Formulation 


Earth-Moon rotating frame Earth-Moon rotating frame 


Mesh Refinement 


0.85 0.9 095 1 1.05 1.1 1.15 


Then solve again with refined mesh 


Traditional SQP Algorithm 


=Minimize the Lagrangian: 
£(#, 4, it) = f@ + 4-hR) + h- G@) 
objective constrain=0O constrain <0 


=This is some nonlinear function which we don’t know how to solve 


"We do know how to solve Quadratic Programming problems, so 
approximate the nonlinear problem as quadratic: 


= Two-term Taylor series expansion of f (X): 
1 
f(%) = f(#*) + F(R") - 6X + 5 OX: Hye" rox 
= One-term Taylor series expansion of constraints: 
h(x) = h(x") + VAR) - 6x 
G(X) = g(x") + Vg(X) - 6x 
=Sequential Quadratic Programming 


= Solve a sequence of quadratic programming (QP) problems that 
approximate the general nonlinear programming problem 


SQP Algorithm Variant 


=Minimize: 
N 3 
f=) duo +540 
t=1 j=1 
Truly quadratic objective 
=Subject to: 
= Dynamics constraints: 
d+]J-5X =0 
d = defects 
ad 
J = 5R 


6X = update to all optimization variables 


Endpoints 


=Previously demonstrated that this approach (or even ordinary least 
squares) can be used to optimize trajectories when the endpoints and 
time of flight are fixed 


=Now extend to variable endpoints & time of flight 


"Easy (fast) to solve problems with: Earth-Moon rotating frame 
= Linear equality constraints 
=" Quadratic inequality constraints 


L, halo 


= Quadratic cost 
orbit 


"Hard (slow) to solve problems with * 
any higher order 1. 
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=Solution: use linear equality constraints 
and add quadratic endpoint term to cost 


Endpoints 


"Define endpoint: g(t) =L, halo orbit, defined by a set of points in a text file 


> i 
"True endpoint constraint: h, = 


e 
"Quadratic expansion of endpoint: 


k “ate 1 0*q(t) 
G(t) = G(t*) + = gt 2 ae 


"Linear endpoint equality constraint: 
he = | - q(t") += 70) 
Ve 
"Add to objective function: 


f= foath ane ag |e 1) 


=With 6 too small, solution bounces 
around optimal Tt indefinitely 
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"With 6 too large, solution converges 
prematurely on sub-optimal T 
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Earth-Moon rotating frame 


L, halo 
orbit 
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Endpoints 


What if we use a linear expansion with a different set of parameters? 


Tried Modified Equinoctial Elements, unsuccessful 


Works well sometimes (when far from singularities) 


Totally fails sometimes (when close to singularities) 
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Endpoints 


What if we use a linear expansion with a different set of parameters? 


Tried Modified Equinoctial Elements, unsuccessful 


Works well sometimes (when far from singularities) 


Totally fails sometimes (when close to singularities) 


Quadratic 
expansion of 
MEE’s 


Line Search 


=Each solution to the QP problem gives us an update 6X to all 
optimization variables 


Xk+1 = XK 4 @. 6x 


"If the problem is sufficiently linear, the QP update is accurate enough to 
assume a= 1 


=Why do a line search? 
=" We do not trust the solution to the linearized problem 


a ¥K+1 — Xk 4G. 6X 


"For short transfers (<1 revolution), no need to perform line search — the 
problem is sufficiently linear to converge quickly with full steps 


Maratos effect 


Line Search 


=A comment on parameterization 


= Line search is only necessary as the solution takes on 
more revolutions 


= With a different parameterization (i.e. orbital elements), 
the revolutions can be “unwound” to keep the problem 
more linear 

= However, the optimization algorithm is too “smart” for 
this 
= Every orbital element set has some singularity (or multiple) 


= Optimization algorithm will exploit the singularity to find a non-physical 
solution with very low cost 


Example applications 


=Now, two examples, with CRTBP dynamics 
= DRO (distant retrograde orbit) to L, halo orbit 
= DRO to different DRO 


"Initial guess is random 


"Endpoints and time of flight are variable, but only allowed 
to change a small amount each iteration, to preserve 
accuracy of linearization 


Example applications 


Iteration 0 


Example applications 


Iteration 2 


Example applications 


Iteration 8 


Example applications 


Iteration 19 
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This transfer requires 15 days and an acceleration of 1.7E-4 m/s? 
(equivalently, 170 mN for a 1000 kg spacecraft) 


Example applications 
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This transfer requires 29 days and an acceleration of 2.8E-4 m/s? 
(equivalently, 280 mN for a 1000 kg spacecraft) 


Example applications 
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This transfer requires 28 days and an acceleration of 2.8E-4 m/s? 
(equivalently, 280 mN for a 1000 kg spacecraft) 


Fuel Optimal Solutions 


=Previously demonstrated we can easily 
transition from one objective to 
another 

= flul? dt 


= With p = 2, large radius of convergence 


= With p = 1, small radius of convergence 


Thrust magnitude (N) 


= Use homotopy method with control law 
to transition from p = 2top=1 


Time (days) 


Example for Earth-Mars low-thrust 
rendezvous 


Implementation notes 


"Implemented in Julia language, with JUMP optimization 
toolbox and Gurobi as QP optimizer 


"Computation time (40-100 nodes): 


= Each iteration: 
= Set up QP problem: 0.2 — 0.5 seconds 
= Solve QP problem: 0.2 — 0.5 seconds 
= Line search: 0.2 — 0.5 seconds 


= Short transfers total time 
=" From random initial guess: 10 — 30 seconds 
=" From close initial guess: ~1 — 3 seconds 
= Long transfers total time varies 
= Line search becomes necessary, so more iterations required 
= Does not always converge 


Low-Thrust Trajectory Optimization 
with Moditied SQP Algorithm 


Questions? 
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